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We show the conclusions claimed in the manuscript arXiv:1202.5309vl by Cuoco, Komatsu and 
Siegal-Gaskins (CKS) are not generally valid. The results in CKS are based on a number of simpli- 
fying assumptions regarding the source population below the detection threshold and the threshold 
flux itself, and do not apply to many physical models of the blazar population. Physical blazar popu- 
lation models that match the measured source counts above the observational threshold can account 
for ~60% of the diffuse gamma-ray background intensity between 1-10 GeV, while the assumptions 
in CKS limit the intensity to <30%. The shortcomings of the model considered in CKS arise from 
an over-simplified blazar source model. A number of the simplifying assumptions are unjustified, in- 
cluding: first, the adoption of an assumed power-law source-count distribution, dN/dS, to arbitrary 
low source fiuxes, which is not exhibited in physical models of the blazar population; and, second, 
the lack of blazar spectral information in calculating the anisotropy of unresolved gamma-ray blazar 
emission. We also show that the calculation of the unresolved blazars' anisotropy is very sensitive to 
the spectral distribution of the unresolved blazars through the adopted source resolution threshold 
value, and must be taken into account in an accurate anisotropy calculation. 



INTRODUCTION 

The contribution of unresolved blazars to the diffuse 
gamma-ray background (DGRB) has been of interest for 
some time (see Ref. [1], hereafter ABH, for a discussion). 
The recent manuscript by Cuoco, Komatsu, and Siegal- 
Gaskins |2] (hereafter CKS) has derived limits on the con- 
tribution of blazars to the DGRB from a combination of 
measurements of the DGRB anisotropy, source-count dis- 
tribution, and intensity. Using a simplistic dN/ dS for the 
blazars, and neglecting any blazar spectral information, 
the CKS analysis concludes that blazars can contribute 
no more than 30% of the DGRB intensity, independent of 
the measured angular correlation power in the DGRB. In 
this note, we show that the CKS limit on the blazar con- 
tribution to the DGRB intensity is not generally valid, 
and strictly the result of their chosen over-simplified 
model. Such a model neglects many crucial features of 
physically- motivated blazar models, e.g. Refs. [3 |3HZ]- 
Importantly, using a physically-constrained source-count 
distribution above the threshold that is consistent with 
that measured by the Fermi-LAT collaboration [8] (here- 
after FBIO) and assumed by CKS, ABH find an intensity 
contribution to the DGRB between 1-10 GeV of approx- 
imately 60%, in direct contradiction to the general claim 
in CKS of a required < 30% contribution. In this note, 
we summarize the reasons for this discrepancy, which re- 
side in a number of invalid assumptions in CKS. 

BLAZAR FLUX SOURCE-COUNT 
DISTRIBUTION FUNCTION 

The blazar source count distribution functions dN/dS 
give the number of blazars expected at a particular flux 




FIG. 1. The best-fit model for the source-count distribution 
function dN/ d-Fioo for the ABH blazar model (solid red line) . 
The data are from FBIO. The power-law fit to dN/dFioo from 
FBIO is shown for comparison (dashed blue line). Note that 
in this plot, _Fioo = /j^q ^^^{dN/dE) dE for comparison to 
the data, instead of S as defined in the text. 

S, defined in this note as 
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S= ds'-^, (1) 
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which is consistent with the flux S defined in CKS, 
though iS* can be defined in different energy bands. The 
Large Area Telescope (LAT) on the Fermi Gamma-Ray 
Space Telescope has measured the dN/dS for blazars to 
be consistent with a broken power-law over the Fermi- 
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LAT sensitivity range in FBIO. CKS assumes that the 
fuh blazar dN/ dS follows this faint-end single power-law, 
down to zero flux. They calculate the diffuse blazar in- 
tensity in the 1 — 10 GeV band as 

f^^ dN 

and they determine the value of St , the "flux sensitivity 
threshold" below which all sources are undetected by the 
Fermi-LAT, by normalizing this intensity integral to the 
measured blazar intensity from FBIO. This definition of 
St neglects the strong spectral dependence of the Fermi- 
LAT point source sensitivity. CKS then calculates the 
Poisson term of the angular power spectrum for unde- 
tected blazars as 

Cp = ^S'dS , (3) 

and compares this value to the measured Cp in the 
DATA: CLEANED sample from Table II of Ref. ^. Using 
the limits on / and Cp, they conclude that "unresolved 
blazars account for only 30% of the IGRB intensity but 
100% of the angular power. However, the CKS analysis 
makes several simplifications which change these results 
greatly when examining physical blazar models. Impor- 
tantly, CKS incorrectly use the same value of St in calcu- 
lating Cp, equation ([s]), as used in /, equation even 
though the point source exclusion limit for the former 
is the IFGL catalog, TS = 25, and that for the latter 
DGRB intensity is more conservative, TS = 50. 

Unlike the broken power-law dN/dS used by CKS, at 
low fluxes dN/dS is expected to flatten rather than con- 
tinuing to increase down to zero flux. The blazar model 
of ABH, for example, exhibits a flattening of source-count 
distribution at low fluxes, as shown in figure [T] The ABH 
model was determined by using a luminosity-dependent 
density evolution (LDDE) model of the gamma-ray lu- 
minosity function using the spectral energy distribution 
(SED) sequence for blazars, which was fit to the Fermi- 
LAT source-count distribution of FBIO and the Fermi- 
LAT-measured DGRB of Ref. PI5] using the spectrally- 
dependent flux limit of sources. Rather than continu- 
ing to rise as a power-law for low fluxes, this source- 
count distribution flattens and provides less blazars at 
low flux than a simple power-law extrapolation would 
predict, as is shown in figure [TJ Additionally, a power- 
law extrapolation down to zero flux of the type assumed 
in CKS is mathematically inconsistent, giving a diver- 
gent number of blazars within our cosmological horizon, 
N — J^^ {dN/dS)dS, while the physical source-count dis- 
tribution of ABH does not. 



Aside from a change in the number of sources at low 
flux, the definition of the flux itself leads to a dN/ dS at 
low fluxes which is different from the broken power-law 
used by CKS. Consider the intensity coming from near 
the threshold flux St: Si < St < 5*2 

f^^ dN 

/band = dS ^'^^ ■ ^^'^ 

Let us consider the sources' fluxes to be from a popu- 
lation that changes its spectrum and/or number density 
from a distribution A^i to N2 near the threshold flux, 

f J^^ {dNi/dE) dE for S < St 
S=\ ^ . (5) 

{dN2/dE) dE for S > St 

The flux from the sources is spectrally dependent, and 
could be, e.g., two different average power law spectral 
distributions 

The intensity calculation therefore changes its form from 
just below to just above the threshold, and the former 
lacks any dependence on the latter, 

/■•^' dNi f^^ dNo 

Therefore the contribution from below the threshold has, 
in general, complete independence from that above the 
threshold, and it should be clear that the contribution 
from just below the threshold can be very different from 
just above the threshold. Therefore, all predictive power 
of equation^ (and therefore in CKS) is lost. The same 
problem can be illustrated with equation [3j which was 
used by CKS to calculate Cp. 

Equation [2] is only valid when blindly extrapolating a 
power-law for the source population to fluxes below the 
observable flux threshold, and furthermore requires the 
assumption of an invariant source population spectrum 
below that threshold. Though neither of these assump- 
tions is valid in general blazar source population mod- 
els, CKS has adopted both. What is properly needed in 
understanding the nature of the source population be- 
low the threshold flux is a physical picture of the blazar 
population. One such physical model is provided by the 
LDDE plus SED sequence model in ABH. 

DIFFUSE INTENSITY, ANGULAR 
CORRELATION, AND THE THRESHOLD FLUX 



^ CKS refers to the DGRB as the isotropic diffuse gamma-ray 
background (IGRB). 



In addition to the simplistic assumptions used in 
extrapolating dN/dS below the Fermi-LAT sensitivity 
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threshold, CKS mishandles the calculation of the thresh- 
old itself. In CKS, the threshold flux is calculated 
using equation [2] and normalizing / to the measured in- 
tensity reported in FBIO. The calculation of Cp is then 
made using the St calculated from /. However, Cp is 
highly sensitive to barely-unresolved sources near the 
threshold, so the calculation of Cp is strongly dependent 
on the chosen value of 5't. A factor of two change in St 
only changes / by ~ 20% but can change Cp by a factor 
of three. 

CKS additionally considers the blazar model of Ref. [5] 
and calculates Cp for this model. However, they use 
the St previously calculated for the FBIO blazar model, 
which was normalized to a significantly different value 
of / than Ref. [S] calculates. This flux threshold value 
is not the correct one for the Ref. [5] blazar model, 
and therefore, the CKS-calculated value of Cp for this 
model is not valid. For comparison, for the ABH blazar 
model, we find that the blazar intensity from 1-10 GeV 
is 2.2 X 10^^ ph cm~^ s~^ sr"-'^, approximately 60% of 
the DGRB. For this intensity, we find a flux threshold of 
St = 2.9 X 10~^ ph cm"'^ s~^, which is much different 
than the = 3.7 x 10-^° ph cm'^ s'^ from CKS. 

Additional problems with CKS are related to equa- 
tions [2] and |3] for the intensity and anisotropy of the 
source population. There is an inherent integration and 
averaging over the source spectrum in these expressions. 
Figure 1 of FBIO shows the threshold flux for the Fermi- 
LAT to be not a single flux value but rather a strong 
function of the blazar spectral index. Depending on the 
blazar index, the threshold flux can vary by two orders of 
magnitude. This is important because the blazar inten- 
sity / is more sensitive to hard sources than the blazar 
anisotropy Cp, so the threshold flux St is, in general, not 
the same for the calculation of / and the calculation of 
Cp. 

As a simple example, the blazar model of FBIO, which 
extrapolates the Fermi-LAT broken power-law dN/dS 
below the Fermi-LAT threshold, considers a spread in 
the blazar spectral indices F: 



dN 
dSdV 

fiS) 



g(r) = exp 



fiS) 9{T) 

A S~^' 

S-P2 



S>Sb 
S<St, 



2a2 



(8) 
(9) 
(10) 



Including the blazar index distribution, equations [2] and [3] 
become 



,00 .s^m 

1=1 dTl dSS — 



-00 Jo 



dN 



(11) 



(12) 



Using the 0.1 — 100 GeV band model from table 4 of 



FBIO and the threshold fluxes in FBIO figure 1, we cal- 
culate / and Cp for this blazar model. Note that this 
analysis was done using Fioo = /^qo McV^^N / dE) dE 
rather than 5, to be consistent with FBIO figure 1. The 
FBIO model gives I — 2A x 10~^ ph cm^^ sr~^ and 
Cp = 3.7x lO^i** (ph cm"2 3-1^2 ^^.-i^ Using equations[2] 
andjSj the equivalent index- independent threshold fluxes 
for each calculation are St{I) ~ 1.7 x 10~^ ph cm"^ s^^ 
and St{Cp) = 3.8 x 10^^ ph cm^^ s^^. As shown above, 
the large difference in St significantly affects the calcula- 
tion of Cp. CKS fails to take this effect into account. 

For a full LDDE plus SED sequence blazar model, 
like ABH, the flux-dependence of the blazar spectrum 
must also be taken into account. An LDDE-based blazar 
model which integrates over blazar luminosity and red- 
shift has been considered, but only with a simple power- 
law energy spectrum, rather than the full blazar SED [6^. 
To do an accurate calculation of the LDDE plus SED 
model from ABH, an extension of the calculations of 
Ref. [5] to include the blazar SED must be done |11) . 



CONCLUSIONS 

As shown in ABH and above, using a source-count dis- 
tribution that is consistent with that measured by the 
Fermi-LAT collaboration (FBIO) and assumed by CKS 
above the threshold, ABH find an intensity contribution 
to the DGRB between 1-10 GeV of 60%, in direct contra- 
diction to the general claim in CKS of a required < 30% 
contribution. The CKS calculation of the Poisson term 
of the angular power spectrum for undetected blazars is 
inadequate. The broken power-law dN/dS they choose 
cannot be accurately extrapolated below the Fermi-LAT 
fiux threshold, and doing so leads to unphysical results. 
CKS also use the incorrect value for the threshold flux 
when calculating Cp and comparing model intensity re- 
sults. Furthermore, the model they consider fails to ac- 
count for blazars' spectral properties, which can affect 
the anisotropy calculation significantly. They assume a 
spectrally-independent threshold fiux for the Fermi-LAT, 
which does not match the actual Fermi-LAT measure- 
ments. For the other model CKS considers, from Ref. [S], 
they use a value of the Fermi-LAT flux threshold which 
does not accurately reflect the threshold flux for that 
model, and therefore this model's exclusion by the Cp is 
questionable. 

Forthcoming work should accurately consider the con- 
sistency between angular correlations in the DGRB and 
its intensity, as contributed by blazars in physically- 
motivated blazar models, and should not rely on unjusti- 
fied extrapolations and the other unqualified assumptions 
present in CKS, as described above. 
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